Terahertz radiation driven chiral edge currents in graphene 
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We observe photocurrents induced in single layer graphene samples by illumination of the graphene 
edges with circularly polarized terahertz radiation at normal incidence. The photocurrent flows along 
the sample edges and forms a vortex. Its winding direction reverses by switching the light helicity 
from left- to right-handed. We demonstrate that the photocurrent stems from the sample edges, 
which reduce the spatial symmetry and result in an asymmetric scattering of carriers driven by the 
radiation electric fleld. The developed theory is in a good agreement with the experiment. We 
show that the edge photocurrents can be applied for determination of the conductivity type and the 
momentum scattering time of the charge carriers in the graphene edge vicinity. 



PACS numbers: 73.50.Pz, 72.80.Vp, 81.05.ue, 78.67.Wj 

The "bulk" transport properties of graphene have been 
studied intensively in recent years and yielded insight 
into the half-integer and fractional quantum Hall effect, 
phase-coherent effects or spin transport on the microm- 
eter scale, to name a few examples [l|, \^. While the 
details of each of those effects depend crucially on the 
linear dispersion relation of graphene and its specific 
material properties, most of the transport phenomena 
have already been studied in other two-dimensional sys- 
tems. Graphene edges, on the other hand, were predicted 
to show insulating or metallic, even magnetic behavior, 
depending on the crystallographic orientation and edge 
chemistry. In scanning tunneling experiments, an en- 
hanced edge density of states was shown [Sl \^ and Ra- 
man scattering experiments provided evidence for the de- 
pendence of scattering mechanisms on the edge orienta- 
tion p, |6| . In transport experiments edge effects are usu- 
ally masked by bulk properties, nonetheless the graphene 
edges are expected to play a crucial role in the electronic 
properties of graphene-based nanoscale devices. 

Here, we present an opto-electronic method to uniquely 
distinguish edge from bulk scattering by exploring edge 
photocurrents in graphene samples illuminated by tera- 
hertz (THz) radiation. For circularly polarized light the 
edge current is observed to form a vortex winding around 
the edges of the square-shaped samples. Its direction re- 
verses upon switching the radiation helicity from left- to 
right-handed. Evidently, the photocurrent is caused by 
the local symmetry breaking at the sample edges result- 
ing in an asymmetric scattering of carriers driven by the 
radiation electric field. It gives rise to a directed electric 
current along the sample boundary in a narrow stripe of 
width comparable to the mean free path. We show that 
the photocurrent measurements provide direct access to 



electron scattering at the graphene edges and allow to 
map the variation of scattering times along the edges. 

We investigated two types of single-layer graphene 
samples: (i) large-area epitaxial graphene prepared by 
high-temperature Si sublimation of 4H and 6H poly- 
types of semi-insulating SiC substrates [7H9| and (ii) 
small area exfoliated graphene flakes [l[ deposited on oxi- 
dized silicon wafers. Below, we report results on epitaxial 
graphene samples (labeled #1-4H, #2-4H, and #3-6H) 
and three samples prepared from exfoliated graphene. 
Hall measurements indicate that the epitaxial samples 
are n-doped (due to charge transfer from SiC [7|) while 
the exfoliated samples are p-doped. The measured carrier 
density lies in the range (2 ^ 7)xl0^^ cm~^, the Fermi 
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FIG. 1: Photocurrent in sample #1-411 as a function of the 
angle ip defining the light polarization. Solid line is a fit to 
Eq. ID) [see also Eq. ((Tj) and discussion]. The inset shows 
the experimental geometry. The ellipses on top illustrate the 
polarization states for various ip. 
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FIG. 2: Photocurrent Ja in sample 7^1-4H as a function of 
the laser spot position. The laser spot is scanned along y 
and the current is picked up from two contact pairs at the 
top (red circles) or bottom (blue full circles) sample edges 
aligned along x (see inset). Dashed lines represent the laser 
beam spatial distribution, which is measured by a pyroelectric 
camera, scaled to the current maximum. The bottom inset 
shows the scan for sample #3-6H. 



energy Ep ranges from 200 to 300 meV and the mobil- 
ity is about lOOOcm^/Vs at room temperature. Ohmic 
contacts were made at samples' edges (see, e.g., inset of 
Fig. [T|). Details on the material growth and characteri- 
zation can be found in [lOj . 

The experiments on edge photocurrcnts are performed 
applying alternating electric THz fields of a high power 
pulsed NH3 laser [llMl3j operating at wavelengths A ~ 
90.5 /xm, 148 /im or 280 /xm (frequencies / = 3.3 THz, 
2 THz and 1.1 THz, respectively). The radiation induces 
indirect (Drudc-likc) optical transitions, because the pho- 
ton energies arc much smaller than the carrier Fermi en- 
ergy. The NH3 laser generates single pulses with a dura- 
tion of about 100 ns, peak power of P w 10 kW, and a 
repetition rate of 1 Hz. A typical spot diameter from 1 
to 3 mm. The beam has an almost Gaussian form, which 
is measured by a pyroelectric camera [IJ] . 

All experiments are performed at normal incidence of 
light and at room temperature. Elliptically and, in par- 
ticular, circularly polarized radiation is obtained apply- 
ing A/4 quartz plates. The resulting pola rization state 
described by the Stokes parameters [l5| Si, S2, and 
S3 = Pcirc is directly related to the angle (p between the 
initial linear polarization of the laser light along the y- 
axis and the plate optical axis. The experimental geome- 
try is shown in Figs. [U [51 and|31 The current is measured 
via the voltage drop across a 50 ft load resistor. 
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FIG. 3: (a) Experimental geometry for the study of edge 
photocurrcnts. (b) Schematic illustration of the edge current 
generation. The electric field of circularly polarized radiation 
rotates clockwise or counterclockwise resulting in a circular 
motion of carriers, which is sketched by red and blue trajec- 
tories, respectively. Our theoretical model, see Eq. 0, shows 
that the circular edge current stems from carriers moving to- 
wards the edge. It is due to the second order .E-field correc- 
tion to the distribution function and involves the retardation 
of the electron motion with respect to the instantaneous elec- 
tric field. Switching the radiation helicity reverses the motion 
direction and, consequently, the electric current, (c) and (d) 
photocurrent topology. Red and blue arrows show the current 
direction for a+ and a- polarizations, respectively. Numbers 
indicate the photocurrent amplitude Ja in microampers. 



Illumination of the edge of unbiased large-area sam- 
ples between any pair of contacts results in a photocur- 
rent. By contrast, if the laser spot is moved toward the 
center the signal vanishes. The detected signal depends 
strongly on the radiation polarization. Fig. [T] The prin- 
cipal observation is that for right- (cr_|_) and left-handed 
((T_) polarizations, i.e., for ip = 45° and 135°, the signs of 
the photocurrent J are opposite. The overall dependence 
J{(p) is more complex. It is well described by 

J{ip) ~ Ja sin 2(y3 + ( Js /2) sin Aip ~ Jc cos^ 2(/3 + ^ 

= JAPcircip) + JBS2{p) + JcSlip)+^ (1) 

and corresponds to the superposition of the Stokes pa- 
rameters with different weights. The first term given by 
the coefficient J a is just proportional to the radiation he- 
licity, whereas the second (J ex Jb) and third (J ex Jc) 
terms change with degree and orientation of the linear 
polarization. Note that the observed offset £, is usually 
smaller or comparable to Ja, Jb, and Jc (see Fig. [T]). 
In our present study we focus on the helicity driven pho- 
tocurrent J A- This is the only contribution which re- 
verses the current direction upon switching the radiation 
helicity from a+ to cr_. We also note that, for circu- 
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FIG. 4: Photocurrent Ja measured at 2 THz for different 
edge segments [see Figs. |3J: and [HJl for circles and triangles, 
respectively]. Lines are fits to Eq. ([7]). The fitting parame- 
ters r/lO"^* s for sample #1-4H and 7^2-4H are indicated by 
numbers. The inset shows the measured circular photocur- 
rent Ja{ojt) at one of the edge segments of sample #1-4H 
(open circles) together with the fit after Eq. ((7]). Data point 
at ojr = 4.4 is obtained applying pulsed CO2 laser. 



larly polarized light (Pcirc = ±1 and 5*1 = 5*2 = 0) and 
^ = 0, the current is solely determined by the first term 
in Eq. ([T]). Therefore, it would be sufficient to measure 
the response to circularly polarized radiation. However, 
to increase the accuracy, we always measured the whole 
polarization dependence, like the one shown in Fig. [TJ 
and extracted Ja by fitting Ec^. Q to the data. 

To prove that the photocurrent is caused by illuminat- 
ing the graphene edges, we scanned the laser spot across 
the sample along the y-axis. The signal was picked up 
from a pair of contacts at the sample top and bottom 
edges aligned along the a;-axis. The experimental geome- 
try and the photocurrent Ja versus the spot position are 
shown in Fig. [31 The current reaches its maximum for 
the laser spot centered at the edge and rapidly decays 
with the spot moving. Comparison of JA{y) with the 
independently recorded laser profile (dotted line) shows 
that the signal just follows the Gaussian intensity pro- 
file. This observation unambiguously demonstrates that 
the photocurrent is caused by illuminating the sample 
edges. Moreover, Fig. [2] reveals that the helicity driven 
current J a changes its sign for opposite edges. 

The above results show that the current direction at 
a specific edge depends on the light helicity. To check 
this in more detail, we investigated currents excited by 
circularly polarized radiation for different pairs of con- 
tacts. Here, the laser spot is always centered between 
the contacts, see Fig. [3^. The current direction for a+ 
(red arrows) and cr_ (blue arrows) circularly polarized 
radiation and the magnitude of J a for various contact 
pairs are shown in Fig. [3b and[3t. The figures document 
a remarkable behavior of the circular edge photocurrent: 
it forms a vortex winding around the edges of the square 



shaped samples which reverses its direction upon switch- 
ing from right- to left-handed. These dependencies are 
observed for all used wavelengths and samples. Helicity 
driven currents have also been observed for small area 
graphene flakes, see [10| for details. 

The observation that a photocurrent occurs only if the 
laser spot is adjusted to an edge agrees with the symme- 
try analysis: In the ideal honeycomb lattice of graphene 
and for our experimental geometry, any photoelectric ef- 
fect is forbidden [16[ , because the two-dimensional struc- 
ture of graphene possesses a center of space inversion. 
Thus, the appearance of photocurrents at normal inci- 
dence of radiation is a clear manifestation of the symme- 
try reduction of the system, in our case, due to the edges. 
We also note that the typical photon energy fiio^lQ meV 
used in experiment is much smaller than the characteris- 
tic energy of carriers Ep ^ 100 meV. Thus, the mecha- 
nism of current formation can be treated classically and 
should involve the action of the light's electric field on 
free carriers in the vicinity of a graphene edge. 

A microscopic process actuating the edge photocurrent 
generation is illustrated in Fig. ^p. It involves the time 
dependent motion of carriers under the action of the elec- 
tric field of circularly polarized radiation and scattering 
at the sample edge. We note that this mechanism is simi- 
lar to that of the surface photogalvanic effect observed in 



bulk materials [17|, [18| . The microscopic theory of edge 



currents is developed in the framework of the Boltzmann 
kinetic equation. In this approach, the electron (hole) 
distribution is described by the function /(p, x, t). It de- 
pends on the carrier momentum p, coordinate x {x > Q 
for a semi- infinite layer), time t, and obeys the equation 
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where E{t) 

radiation, v ~ vp/p is the electron velocity, w « 10^ m/s 
is the effective speed, q is the carrier charge [q = +|e| for 
holes and — |e| for electrons), and Q{f} is the collision 
integral. The distribution function can be expanded in 
series of powers of the electric field, 



/(p, X, t) = foiep) + [flip, x)e' 
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(3) 

where fo{Sp) is the equilibrium distribution function with 
Sp = vp being the electron energy, /i oc \E\, and 
/2 OC \E\'^. The first order in E correction to the dis- 
tribution function oscillates with frequency w and does 
not contribute to a dc current. The directed electric cur- 
rent along the structure edge is, therefore, determined by 
the second order E-iield correction /2 and given by 
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(4) 



The factor 4 accounts for the spin and valley degeneracy. 



We solve Eq. ([2]) and calculate the current, Eq. Q, for 
the simple form of the collision integral, 

Q{/(P,.,t)}=.-fe^i^)-Mp), (5) 

with T being the scattering time, and the boundary con- 
dition at a; = 0, 

fip^ >0,Py,0,t) = 

- I ^fip', 0, msp - Sp.)ei-v',)dp' , (6) 

corresponding to diffusive scattering. In the case of a 
degenerate gas, the edge current takes the form (see [10[ 
for details): 
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The helicity-driven current is given by the first term be- 
cause \[Eq X E^]z = — -Pcirc for our geometry where the 
light propagates along —z. The second term yields the 
current caused by linearly polarized radiation and van- 
ishes for circular polarization. In the case of elliptically 
polarized light, Eo,xE^y + Eo^yE^.^ oc (1/2) siniip = 82- 
Both contributions are clearly detected in the experiment 
and correspond to the first (ex Ja) and second (oc Jb) 
terms in the empirical Eq. ([1]), see Fig. [TJ 

The helicity driven photocurrent described by Eq. ([7]) 
vanishes for zero frequency, has a maximum at ujt ~ 0.6 
and decreases rapidly at higher frequencies. Exactly this 
behavior is found in experiment (see inset of Fig. U) as 
we explain in more detail below. The only free parameter 
in Eq. (O is the scattering time r. Corresponding data 
are shown in Fig. |3] where the photocurrent values mea- 
sured at 2 THz for each of the contact pairs are plotted. 
These data points are first compared to calculated traces 
of Ja employing Eq. ([7|) . Solid lines are calculated using 
the bulk values for the time r extracted from resistivity 
and carrier density for samples #1-4H (r = 2.0 x lO"^** s) 
and #2-4H (r == 2.8 x 10"" s). The bulk scattering times 
used in Eq. ([7]) give for some of the contact pairs already 
perfect quantitative agreement. For other edge segments 
the current deviates significantly. This is a consequence 
of the strongly non-linear dependence of Ja on r. Vary- 
ing r by only ±15% changes the current by ±50 %. By 
fitting the photocurrent Ja we can extract the local scat- 
tering time T for every edge segment shown in Figs. [3J; 
and d). The best fits are shown by dashed lines in Fig.|4] 
and constitute a map of scattering times along the edge. 
The average value of the circular edge current scales with 
the sample mobility. To check the frequency dependence 
predicted by Eq. ([7]) we show in the inset of Fig.U] J^ vs. 
LOT for one edge segment using the extracted r. The data 



points are perfectly described by Eq. ([7]) and confirm the 
model. 

While the magnitude of the circular edge photocurrent 
agrees well with theory, the expected polarity of Ja for 
n-type graphene is opposite to the one observed. This, at 
first glance, surprising result agrees with results from spa- 
tially resolved Raman measurements demonstrating that 
edges of n-type graphene layers exhibit p-type conduc- 
tivity [ay. This explains the sign of the photocurrent, 
which is generated in a narrow edge channel comparable 
to the mean free path (w 10 -^ 20 nm) and has oppo- 
site sign for electrons and holes, see Eq. ^. Actually, 
the difference in the conductivity type can be also un- 
derstood from the details of the sample fabrication. It is 
well established that epitaxial graphene on SiC(OOOl) is 
n-doped due to charge transfer from the interfacial buffer 
layer (see, e.g., [7|,|8[), while so-called quasi- free-standing 
graphene, lacking such buffer layer and sitting on a hy- 
drogen terminated SiC(OOOl) surface, is p-doped jl9{ . 
Therefore, it is reasonable that the edges of epitaxial 
graphene, exposed to the SiC substrate without the in- 
terfacial layer, can be p-doped. This assumption is cor- 
roborated by similar reports on the transition from n- 
to p-type of doping at the edges of graphene flakes on 
Si02, which were attributed to the difference in the work 
functions of graphene and the substrate [20| . 

To summarize, our observations clearly demonstrate 
that illuminating monolayer graphene edges with polar- 
ized terahertz radiation at normal incidence results in 
a directed electric edge current. The effect is directly 
coupled to electron scattering at the graphene edge and 
vanishes in bulk graphene. Our results suggest that cir- 
cular the photocurrents can be effectively used to study 
edge transport in graphene even at room temperature. 
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SUPPLEMENTAL MATERIAL 
SI. Details of the Samples 

We investigated three epitaxial samples grown on SiC. 
Samples #1-4H and #2-4H were grown by the Linkoping 
group on a Si-terminated surface of a 4H-SiC(0001) semi- 
insulating substrate (Cree Inc.) 3- The reaction kinetics 
on the Si-terminated surface is slower than on the C-face 
because of the higher surface energy, which fosters ho- 
mogeneous and well controlled graphene formation [Tj. 
Graphene was grown at a temperature of 2000°C and 
1 atm Ar gas pressure resulting in monolayers of graphene 
atomically uniform over more than 1000 /zm^, as shown 
by low-energy electron microscopy [2l|. Eight contacts 
were produced by depositing Snm of Ti and 100 nm of 
Au. The quadratic sample size of 5 x 5 mm^ was achieved 
by oxygen plasma etching of all four edges. Hall mea- 
surements indicate that the large area samples are n- 
doped due to charge transfer from SiC 0, i, [22l-ii| . The 
measured carrier concentration is between 3x10^^ cm~^ 
and 7x10^^ cm~^, the Fermi energy Ep ranges from 200 
to 300 meV and the mobility is about lOOOcm^/Vs at 
room temperature. In these samples, as well as in other 
large-area samples, the resistance at room temperature 
is about 2 to 5 kil. 

The third epitaxial graphene sample #3-611, was 
grown by the Erlangen group on 6II-SiC(0001) wafers 
(II- VI Inc.). Graphene growth was performed using sub- 
limation growth in Ar atmosphere [7|, |9[. First, polishing 
damage was removed by etching the substrate in 1 bar 
hydrogen at 1550°C for 15 min. Second, graphene was 
grown by annealing the sample in 1 bar Ar at a temper- 
ature of 1650°C for 15 min. The graphene coverage was 
determined by x-ray photoelectron spectroscopy (XPS). 
The square-shaped sample size of 4 x 4 mm^ was achieved 
by mechanical cutting the edges. Both, carrier density 
and mobility in the n-type sample #3-611 are very simi- 
lar to those of the Linkoping samples. 

For all epitaxial graphene samples, low-tcmpcrature 
quantum Hall measurements reveal the high quality and 
homogeneity. 

The exfoliated graphene samples (small-area graphene 
samples 4, 5, and 6) has been prepared from natural 
graphite using the mechanical exfoliation technique [l| on 
an oxidized silicon wafer. The oxide thickness of 300 nm 
allowed to locate graphene flakes in an optical micro- 
scope and to assess their thickness. We checked the re- 
liability of this method using Raman spectroscopy and 
low-temperature quantum Hall measurements on similar 
samples [26|. The single layer graphene flakes obtained 
by this method were typically p-doped by adsorbed con- 
taminants with carrier concentrations p < 2 x 10^^ cm~^. 
The Fermi energies were Ep < 165 meV and the mobili- 
ties at room temperature of the order of 2.5x10^ cm-^/Vs. 



The flakes included in this study were all single layer with 
the flakes size of the order of 10 to 30 micrometers. The 
sample morphology was characterized by atomic force mi- 
croscopy measurements under ambient conditions with 
the microscope in intermittent contact mode with stan- 
dard silicon tips [27| . After recording the position of the 
flakes with respect to predefined markers, we contacted 
them by electron beam lithography and thermal evapora- 
tion of 60 nm Pd electrodes. The resistance of graphene 
measured between various contacts was about 1 to 3 kJ7. 



S2. Laser beam parameters 

In addition to the pulsed THz laser described in the 
main text, we also used a continuous-wave (cw) CH3OH 
laser (A = 118 /^m) with a power of P w 20 mW. In the 
experiments applying the CH3OH laser, the cw radiation 
was modulated at chopper frequencies in the range from 
120 to 600 Hz. The sign of the signal is defined as a 
relative phase with respect to the lock-in reference signal, 
which was kept the same for all measurements. 

Elliptically and, in particular, circularly polarized radi- 
ation has been obtained by transmitting the laser beam, 
which is initially linearly polarized along the y-dircction 
for the pulsed laser and along the x-dircction for the cw 
laser, through A/4 crystal quartz plates. The resulting 
polarization state is directly related to the angle f be- 
tween the initial linear polarization of the laser light and 
the optical axis of the plate. It is described by the Stokes 
parameters 6*1, ^2, 5*3 [la]. In particular, the dependence 
of the circular polarization degree, given by S3, on the 
angle ip in our experimental geometry has the form 



Ph 



5*3 ((p) — sin2iy9 , 



(8) 



The parameters Si and 5*2 are given by the bilinear com- 
binations of the polarization vector components. 
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The parameters 6*1 and S2 describe the degree of lin- 
ear polarization in the coordinate axes a;, y and in the 
coordinate frame rotated about an angle of 45° , respec- 
tively 15|. Note that radiation is incident along —z axis. 
The resulting polarization ellipses for the cw THz laser 
for some angles (f are sketched on top of Fig. [5l Further 
details on the experimental technique can be found in 
Ref. ii. 



S3. Photocurrents in small-area samples 

Helicity driven photocurrents excited at normal inci- 
dence have also been observed in small-area graphene 
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FIG. 5: Photocurrent in small area sample #4 as a function 
of the angle tp. Squares and circles show the photocurrent 
picked up from the contact pairs at the upper (Ju) and lower 
(Ji) edges of the graphene flake, respectively (see inset). Full 
lines show fits to the calculated total current after Eq. (|10|l . 
The inset demonstrates the experimental geometry, Uu and 
Ui are voltages recorded by lock-in amplifier. Arrows show 
current directions for right circularly polarized light. The 
signal is excited by radiation with the wavelength A — 118 ^m, 
power P ~ 20 mW and a diameter of the laser spot about 
1.5 mm. On top, the polarization ellipses corresponding to 
various angles <p are illustrated. 



flakes. Examples of the current polarization dependence 
are shown in Fig. [S] for two different pairs of contacts. 
Similar to data obtained in large-area samples, it can be 
well fitted by Eq. (1) of the main text, which reads 

Ji^) = JAPcn-civ) + JbS2{^) + JcSliip) + e . (10) 



While photocurrents are observed in both large-area 
and small-area samples, the analysis of the edge pho- 
tocurrents is much easier in the large-area samples. Ac- 
tually only in the latter type of samples the illumination 
of a single edge by THz radiation could be realized in 
our experiments. Such selective excitation has enabled 
the accurate analyses of the edge currents. By contrast, 
in micron-sized exfoliated samples the spot size is much 
larger than the graphene fiakes and the effects of different 
edges are superimposed. 



S4. Theory 

Equation (7) of the main text is obtained by expanding 
the distribution function f{p, x, t) in series of the electric 
field. To first order in the electric field, solution of Eq. (2) 



with the boundary condition (6) has the form 
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where /q = dfo{e)/de. The equation for the second-order 
correction /2 (p, x) is given by 
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which yields 
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By using Eq. (4) for the edge electric current and Eq. P^ 
we derive 
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where the above two contributions to the current stem 
from the first and second terms on the right-hand side of 
Eq. P^ . respectively. Finally, taking into account that 
the electron energy and velocity in graphene are given by 
Ep = v\p\ and V ~ vp/p, respectively, and assuming that 
the electron gas is degenerate and r is independent of ep, 
we obtain Eq. (7) of the main text. 

The current Eq. (7) of the main text is consistent with 
the point-group symmetry Cg containing the mirror plane 
{x,z). It should be noted that for elliptical polariza- 
tion, in some experiments the photocurrent described by 
Eq. (7) is superimposed with an additional contribution 
proportional to cos^ 2ip [see the third term in Eq. ([T0|) ]. 
This term can be attributed to a lowering of the system 
symmetry to Ci showing the non-equivalence of y and 
~y directions, e.g., due to (i) inhomogeneous photoex- 
citation, (ii) macroscopic roughness of the investigated 
edges, (iii) non-equivalence of contacts, etc. The sym- 
metry reduction hinders the edge photogalvanic currents 
under study and complicates their analysis. However, the 
obstacle can be easily overcome applying circularly po- 
larized radiation, like used in the present work. Indeed, 
as addressed above the edge currents driven by circularly 
polarized light change their signs upon variation of the ra- 
diation helicity. By contrast, other current contributions 
caused by the additional symmetry lowering are insen- 
sitive to the radiation helicity and vanish for circularly 
polarized radiation. 
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